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The final-state spectrum of f3 decaying tritium anions T ~ was calculated. The wavefunctions de- 
scribing the initial T ~ ground state and the final 3 He states were obtained by the full configuration- 
interaction method. The transition probability was calculated within the sudden approximation. 
The transition probability into the electronic continuum is extracted from the complex-scaled resol- 
vent and is shown to converge for very high-energies to an approximate analytical model probability 
distribution. 
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I. INTRODUCTION 

The neutrino rest mass is a very important parameter 
for cosmology, astrophysics, and the standard model of 
elementary particles. The existence of neutrinos, already 
postulated by Pauli and put into a mathematical frame- 
work of /3 decay by Fermi [l[ long time ago, was verified 
by Reines and Cowan in 1956 [2]. However, despite the 
high solar neutrinos flux of about billions per m _2 s _1 on 
earth the answer to the question about their rest mass is 
one of the big unknowns in physics. Since neutrino-flavor 
oscillations have been observed in the late nineties at the 
Super-Kamiokande experiment 0] a non-vanishing neu- 
trino rest mass has to be expected. Unfortunately, this 
type of experiments reveals only mass differences between 
neutrino flavors. 

The presently constructed KATRIN (Karlsruhe tri- 
tium neutrino-mass) experiment with an expected sensi- 
tivity of about 0.2 eV (90% C.L.) should have the ability 
to determine the absolute value for one of the flavors or 
at least a new upper limit to it [4|. This so called next- 
generation tritium /3-decay experiment is only based on 
kinematic relations and energy and momentum conser- 
vation. Thus KATRIN provides a model-independent 
direct measurement of the antineutrino rest mass mp e 
(more accurately the mass of the antineutrino in a given 
mass-flavor mixture mostly attributed to the electronic 
neutrino). In more detail, will be extracted in a fit 
procedure from the shape of the [3 spectrum. Besides the 
precise measurement of the /3-electron energy spectrum 
it is crucial for the mass extraction to know how the /3 
spectrum is modified by the final-state spectrum of the 
decay product. As in the previous most recent tritium 
neutrino-mass experiments in Mainz and Troitsk the T2 
molecule is chosen as tritium source. T 2 comprises a 
compromise between experimental accessibility and the- 
oretical treatability. The final-state spectrum of its decay 
product 3 HeT + was therefore subject of a number of very 
detailed calculations finally accumulating in the 

one covering the whole energy regime [l4j . Recently, the 
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spectrum was further adapted to specific needs (isotope 
distribution and temperature) of the KATRIN experi- 
ment fisl ]. 

Although a high purity of the molecular tritium source 
is expected for KATRIN, the produced /3 electrons can in- 
teract with other gas molecules and thus produce tritium 
species different from T2 . One of the expected processes 
is the dissociative attachment 

e" +T 2 -> T~ +T (1) 

where T ~~ formation occurs. Despite the relative small 
cross section compared to, e.g., the one for vibrational 
excitation of the T 2 molecule, this process is very impor- 
tant. The reason is the higher endpoint energy of the /3 
spectrum for the decay of T ~~ compared to the one of 
T 2 . Due to this fact the occurrence of T ~ ions leads to 
a systematical error and hence to a possible limitation of 
the sensitivity of KATRIN, if it is not properly accounted 
for 0. 

To the authors' knowledge, there exist so far only two 
theoretical predictions for the final-state spectrum fol- 
lowing ft decay of T _ . However, in [Tfl ] only transition 
probabilities to 4 final states are reported. Furthermore, 
the results in [l7j disagree substantially from the ones 
given in an earlier work [l~8j that was, however, also lim- 
ited to 10 final states. The aim of this work was thus 
to provide a complete final-state spectrum for the decay 
process 

T - ->■ 3 Hc + c" + V c . (2) 
and to shed some light on the disagreeing earlier results. 

II. METHOD AND COMPUTATIONAL 
DETAILS 

The calculation of the non-relativistic eigenstates of 
the atomic systems is performed within the approxi- 
mation of an infinitely heavy mass of the nuclei, i. e. 
T « °°H and 3 He « °°Hc . This is justifiable due to 
the large mass difference of the nucleus and the elec- 
trons. The calculation of the final-state spectrum can be 
performed analytically for neutral T atoms and reveals a 



2 



negligible mass dependence. Hence a large mass depen- 
dence is also not expected in the case of tritium anions. 
The non-relativistic Hamiltonian for the two-electron sys- 
tem has the form (atomic units with m e = l,e = l,ft=l 
are used throughout, if not specified otherwise): 



H = ~(Ai + A 2 )-Z(- + i 
2 V r i r 2 



1 



(3) 



where Z is the charge of the nucleus and the posi- 
tion vector of the i th electron. Due to the fact that the 
only bound state of T _ is a singlet state with angular 
momentum L = AI = [T^ ] , only symmetric spatial con- 
figuration state functions (CSF) are important, 
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To determine the eigenstates and corresponding energy 
eigenvalues a simple expansion in Slater-type orbitals 
(STO) is used, 
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The n, I, m are integer parameters with limitations anal- 
ogously to the ones for the hydrogen quantum numbers 
and the Y™ represent the spherical harmonics. The Q are 
positive real parameters. An appropriate choice of these 
parameters allows for the achievement of an in principle 
complete coverage of the Hilbert space of the one-particle 
part of Hamiltonian (J3]). In the full configuration- 
interaction (CI) method the eigenstates are expressed as 
a linear superposition of all possible symmetry-adapted 
CSFs 
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that can be formed with the aid of the chosen STO basis. 
The expansion coefficients Cjk are determined by solving 
the generalized eigenvalue problem obtained from insert- 
ing the wavefunction ansatz of Eq. (j6]) into the eigenvalue 
equation of the Hamiltonian ([3]). 

The final-state spectrum of He is calculated within 
the sudden approximation [20j that is based on the fact 
that the escaping /? electron has a much higher veloc- 
ity than the bound electrons. In the analysis of tritium 
neutrino-mass experiments like KATRIN only the (3 elec- 
trons with an energy near the endpoint of the j3 spectrum 
at 18.6 keV are used. Their velocity is clearly much larger 
than the average speed of the bound electrons in T ~~ . In 
fact, the validity of the sudden approximation has been 
demonstrated for T 2 in [l(| HH, HH where the first-order 
correction terms were derived and explicitly calculated. 
From those results it is apparent that also for T ~ the 
sudden approximation is expected to be valid within the 
accuracy required for the analysis of an experiment like 



KATRIN. Nevertheless a brief discussion of possible ef- 
fects beyond the sudden approximation on the final-state 
spectrum is given at the end of this work. 

A basis set of 555 STOs yielding 3481 CSFs in the 
full CI calculation was used to obtain the final results 
shown in this work. This STO basis set contains all pos- 
sible kinds of orbitals (with restrictions on I and m as 
mentioned above) up to the angular quantum number 
I = 7, —7 < m < 7. For the optimization of the pa- 
rameters Q a genetic and several other algorithms [22j 
were tested . However, none of those algorithms lead to 
completely convincing results. Therefore, the parame- 
ters were finally optimized by hand. The difficulty of the 
parameter optimization is due to the requirement to con- 
struct a basis set with high coverage of the Hilbert space 
while avoiding inaccuracies due to numerically caused lin- 
ear dependencies. With the aim to achieve a uniform de- 
scription of the possible He final states it is favorable to 
obtain a homogeneous and a high density of states in the 
continuum as well as a large number of bound states. If 
a large number of CFSs is used, the optimization of the 
individual Q values becomes less important, since the 
full CI method leads to a sufficient mixing of the Hilbert 
space covered by the various STOs. Therefore, the pa- 
rameters Q were chosen to start in an interval between 2 
and 3 and to decrease in value for increasing n (for a given 
I). This procedure avoids numerical problems and allows 
the construction of a huge, but linearly independent ba- 
sis set. This basis set is used for both the ground state of 
T ~ and all final states of He . The chosen basis leads for 
T~ to the ground-state energy = -14.3602 eV that 
is only 0.8 meV above the very accurate values in |23[ |24| . 
In the case of He the adopted basis set yields 16 states 
below the ionization continuum. Out of those 16 states 
15 are identified as true physical states, while the 16th 
state is a pseudo state that resembles the remaining in- 
finite number of Rydberg states as a consequence that a 
finite basis set is adopted. 

Within the sudden approximation the transition prob- 
ability for T ~~ decays into bound states of He is simply 
given by the squared overlap 
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of the initial state ), i. c. the T~ ground state, and 
the final state l^^ ), i- e. the n-th bound state of He . 

To calculate the transition-probability density into 
continuum states the complex scaling method is used. 
It is based on the mathematical development by Aguilar, 
Balslev, and Combes [26l . |27| as well as Simon [28| The 
application of this method leads in practice to a sim- 
ple but powerful modification of the Hamiltonian H in 
Eq. ©, 



H(0) = exp(-2i0)T + exp(-i0)V . 



(8) 



In Eq. ([5]) T and V are the usual kinetic and potential en- 
ergy operators of He, respectively. The complex-scaling 
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TABLE I: Population probabilities P n of the S bound states 
of helium after the j3 decay of a T ~ anion. Also given are the 
corresponding energies E n (in atomic units) obtained in the 
present work. 



n 


E n 


Pn(%) 


Pn(%) 


Pn(%) 




(this work)" 


(this work) 


(Ref. IT) 


(Ref. 18) 


1 


-2.9034572 


22.98998 


22.993764 


19.147 


2 


-2.1459527 


46.86960 


46.867404 


21.149 


3 


-2.0612659 


0.01320 


0.135 


0.27 


4 


2.0335841 


0.18363 


0.21 


0.143 


5 


-2.0211749 


0.09220 


— 


0.07 


6 


-2. 0145604 


0.05262 


- 


0.039 


7 


-2.0106235 


0.03275 


— 


0.024 


Q 
O 




r\ noi 
U.Uzl / 




U.UlD 


9 


2.0063679 


0.01522 




0.011 


10 


-2.0051407 


0.01111 




0.008 


11 


-2.0042355 


0.00848 






12 


2. 0035507 


0.00666 






13 


-2. 0030205 


0.00513 






14 


-2.0025951 


0.00469 






15 


-2.0022570 


0.00273 










70.30975 


70.206168 


40.877 


"The bold digits agree 


with the results 


in Ref.|H 





angle 9 can in principle be chosen arbitrarily within 
0° < < 45°. In the limit of an infinite basis all observ- 
ables calculated with the aid of complex scaling should 
become independent of 9. Since only finite basis sets can 
be applied in practice, only approximate eigenstates can 
be obtained that may depend on 9. The angle 9 can thus 
be understood as a variational parameter that modifies 
the adopted basis as can be seen from the inverse rela- 
tion between basis-set exponents and the scaling angle 
discussed, e.g., in (29|. A diagonalization of the Hamil- 
tonian (jHJ) in the basis described by the Eqs. (H|) and §5$) 
yields the complex-scaled energies Ej(9) and wavefunc- 
tions ^j(9) where the latter are still defined by Eq. ©, 
but with complex coefficients Cjk(9). 

With the aid of the complex-scaled energies and wave- 
functions the transition-probability density into the elec- 
tronic continuum can be extracted from the complex- 
scaled resolvent according to [1] 



P(E,9) = 
1 



-Im ■ 



E 



(9) 



The (^>(9*) \ is the biorthonormal eigenstate to \^(9)}. It 
is obtained from the latter by a transposition and com- 
plex conjugation of the angular part, while the radial part 
is only transposed but not complex conjugated. The sum 
over k includes all complex-scaled eigenstates and eigen- 
values calculated by solving the generalized complex sym- 
metric, but non-hermitian eigenvalue problem. As dis- 



cussed above in the limit of exact eigenstates the density 
P(E, 9) becomes independent of the complex-scaling pa- 
rameter 9. A variation of 9 for approximate eigenstates 
provides the possibility to determine an optimal 9 op t with 
highest stability. The best approximation of P(E, 9) is 
then obtained according to 



dP(E,( 
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min. -+ P(E) := P(E,6 opt ). (10) 



Furthermore, the 9 dependence of the spectra gives an 
indication for the convergence of the results. 



III. RESULTS 

In Table U the calculated transition probabilities for 
15 1 S bound states of He are listed. The results reveal 
that almost every second T ~~ decay will end in the first 
excited state of He . The next probable final state is 
the He ground state with nearly 23%. With a summed 
probability of 0.45% the higher excited He states are 
rarely populated after (3 decay of T — . The sum over 
all calculated bound states yields 70.3%. The summa- 
tion over all calculated states (discrete and discretized 
continuum states) yields the expected value of 100.00%, 
since the same basis is used for initial and all final states, 
but indicates the proper numerical implementation. The 
excellent agreement of the energy eigenvalues at the or- 
der of /xhartree with the very accurate data in [25| as- 
sures on the other hand the high quality of the basis set 
adopted in the present work and its ability to describe 
many states simultaneously with high precision. A closer 
view on the energies shows that the degree of accuracy of 
the present results follows the expected trends. First, the 
accuracy increases with n, since the importance of cor- 
relation decreases, if the state becomes more asymmet- 
ric and the two electrons have smaller spatial overlap. 
For even higher values of n the states become increas- 
ingly diffuse and thus it is very difficult to describe them 
properly without running into numerically caused linear 
dependencies. 

A comparison to the final-state probabilities reported 
by Frolov [TtJ and Harston and Pyper [l8[ is also given 
in Table Q] Especially for the highly populated ground 
and first excited states the results of this work confirm 
the expectedly very accurate results of Frolov [17| that 
were obtained with explicitly correlated basis functions. 
The agreement for the third excited state (n — 4) is, 
however, less good, and for the second excited state (n = 
3) there is even an order of magnitude difference. All 
attempts to improve the basis set for this state failed to 
yield a better agreement. This could be an indication for 
a typographical error (a missing zero after the decimal 
point) in |17| . 

The comparison with the results in [181 ] that were 
obtained with a relativistic MCDF (multi-configuration 
Dirac-Fock) method shows on the other hand pronounced 
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TABLE II: Discretized final-state probability distribution P(Ei) for He following the ft decay of a T anion. The mean 
excitation energies Ei are given relative to the ground state of 3 He . 



Ei(eV) 


P(Ei)(%) 


Ei{eV) 


P(Ei)(%) 


Ei(eV) 


P{Ei){%) 


Ei(eV) 


P(Ei){%) 


25.084 


0.36869 


50.096 


0.07054 


75.086 


0.09926 


186.05 


0.00960 


26.081 


0.32908 


51.097 


0.07403 


76.085 


0.09172 


191.06 


0.00873 


27.081 


0.28924 


52.100 


0.08005 


77.085 


0.08480 


196.06 


0.00796 


28.081 


0.25425 


53.104 


0.09041 


78.086 


0.07863 


201.06 


0.00727 


29.082 


0.22462 


54.112 


0.10943 


80.950 


0.31969 


206.07 


0.00666 


30.082 


0.19976 


55.125 


0.14982 


85.967 


0.23208 


211.08 


0.00612 


31.083 


0.17887 


56.158 


0.26863 


90.980 


0.17471 


232.06 


0.03471 


32.084 


0.16126 


57.305 


1.44117 


95.991 


0.13529 


272.24 


0.01998 


33.084 


0.14632 


57.863 


16.80345 


101.00 


0.10714 


313.03 


0.01250 


34.085 


0.13358 


58.911 


0.15914 


106.01 


0.08644 


352.94 


0.00832 


35.085 


0.12267 


59.972 


0.02281 


111.01 


0.07084 


392.90 


0.00581 


36.085 


0.11328 


61.109 


0.00924 


116.02 


0.05883 


433.04 


0.00421 


37.086 


0.10517 


62.095 


3.16161 


121.02 


0.04942 


472.98 


0.00314 


38.086 


0.09816 


63.026 


0.13662 


126.03 


0.04194 


513.67 


0.00239 


39.087 


0.09209 


64.113 


0.13095 


131.03 


0.03590 


553.22 


0.00186 


40.087 


0.08684 


65.091 


0.09512 


136.03 


0.03098 


591.30 


0.00148 


41.088 


0.08233 


66.100 


0.11217 


141.03 


0.02692 


630.58 


0.00119 


42.088 


0.07847 


67.099 


0.12245 


146.04 


0.02355 


654.79 


0.00101 


43.089 


0.07523 


68.105 


0.13639 


151.04 


0.02071 


664.12 


0.00089 


A A AOA 

44.089 


0.07256 


69.171 


0.26130 


156.04 


0.01832 


O f\C\ 

663.99 


0.00081 


45.089 


0.07046 


70.210 


0.20446 


161.04 


0.01629 


666.03 


0.00073 


46.085 


0.06892 


71.109 


0.10126 


166.05 


0.01454 


688.91 


0.00063 


47.084 


0.06803 


72.101 


0.12307 


171.05 


0.01303 


776.21 


0.00048 


48.090 


0.06792 


73.077 


0.12211 


176.05 


0.01173 


1550.1 a 


0.00358° 


49.098 


0.06870 


74.083 


0.10784 


181.05 


0.01059 




29.65399 


"For enerf 


;ies above 904 eV the model tail in 


Eq. (1 1 3lt was used. 











differences. The deviation is most remarkably for the 
first excited state that according to the present work and 
[i~7j should be populated with about 47 % probability and 
thus should clearly dominate the final-state distribution. 
However, in the MCDF results in [l8| its probability is 
found to be about 21%. For the other states, except 
ri = 3, the results in are always smaller than the 
present ones. The deviation increases rather uniformly 
from about 17 to 28 % for n varying between 1 and 10. 
Since relativistic effects are expected to be small for light 
nuclei like T ~ and 3 He , it appears very likely that the 
main reason for the difference of the results in [l8| to 
the present ones (as well as the ones in 17]) is due to 
the small number of configurations used in the MCDF 
method compared with the present full CI method. Un- 
fortunately, no details (like energies) of the MCDF cal- 
culation in [l]l] are available to further clarify this issue, 
but any realistic estimate of the size of relativistic effects 
excludes their responsibility for the large discrepancy be- 
tween the results in [18| compared to the non-relativistic 
calculations of this work or the one in [l7j . 

The calculated transition-probability density into the 
electronic continuum of 3 He is presented for three dif- 




Energy [eV] 

FIG. 1: (Color online) Final-state continuum probability den- 
sity of He after /3 decay of T" for 6 = 24° (red), 30° (blue), 
and 36° (black). (The energy scale is given relative to the He 
ground state.) 
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ferent complex-scaling angles (9 — 24°, 30°, and 36°) in 
Fig. [T] The overall spectrum is practically independent 
of 9. This indicates the high quality of the adopted ba- 
sis set also for describing the electronic continuum. As 
is usually the case, (higher lying) resonances are most 
sensitive to the choice of 9. This is due to the fact that 
it is difficult to find a single value of 9 that is equally 
appropriate for describing a certain resonance and the 
underlying background continuum. 

The continuum probability density is dominated by 
a peak corresponding to the first doubly excited singlet 
state 2s 2 . About 19% of the T~ decays ends up in the 
energy interval between 54.5 eV and 60 eV. Above the 
65.4 eV threshold the higher-lying doubly excited states 
2sns and in the regime up to 79 eV (with diminishing im- 
portance) the 3sns peaks can be identified. The complex- 
scaling method provides the probability density P(E) at 
any value of E and thus as a continuous function. In 
view of the sharp resonant structures and in accordance 
with the experimental needs, the final-state distribution 
is given in a discretized form as in [141 ]. For this pur- 
pose, the probability distribution P(E) has been divided 
into small bins covering an energy range of 1.0 eV (up 
to a transition energy of 78.59 eV), 5.0 eV (from 78.59 to 
214 eV), and 40.0eV (from 214 to 904eV). For each bin 
the average excitation energy E^ and the integrated tran- 
sition probability P{E%) were calculated and are given in 
Table HI 

For the high-energy continuum states (above 904 eV) 
an approximate model tail is introduced, similar to the 
case of T2 [HI, 03 ■ However, the situation is more com- 
plicated for T _ . In T2 j3 decay the high-energy tail was 
derived based on the idea that for sufficiently large ener- 
gies of the escaping (formerly bound) electron the effec- 
tive potential of the remaining 3 HcT 2+ ion can be well 
approximated by a point charge Z = 2. In fact, the 
remaining electron and tritium nucleus may be viewed 
as pure spectators and thus the transition probability 
should approach for high energies the one obtained for a 
/3-decaying tritium atom for which an analytical result is 
known. Due to the existence of two equivalent electrons, 
the atomic result is simply multiplied by a factor of two 

While for T2 a hydrogenic wavefunction is a reason- 
able first-order approximation for the initial state, this 
is not the case for T _ . In fact, within independent- 
particle models T~ is unstable. As a consequence, the 
fast electron in the final state may be well represented 
by a Coulomb wavefunction for a point charge Z = 1 
(formed by the remaining He + ion), but the modeling 
of the initial state is less obvious within an independent 
particle model. This is also evident from the alternative 
point of view that a description of the remaining T nu- 
cleus and bound electron as a spectator would correspond 
for the active electron to an initial state with Z = and 
thus no bound state. In order to obtain an atomic-like 
high-energy tail the initial state is thus modeled as a 
hydrogen-like state with variable exponent. This expo- 



nent is then obtained by fitting the model spectrum to 
the ab initio spectrum of the full two-electron calculation 
in the energy range between 500 eV and 10,000 eV. 

The initial T ~ ground-state wavefunction (omitting 
the spin part for better readability) is then approximated 



as 



i*r> = ii« z,i « 



(ii) 



where (r|ls z =) is an STO with (n,l,m) = (1,0,0), i.e. 
an atomic hydrogen Is orbital with effective charge Z e . 
In the spirit of the sudden approximation the spectator 
electron remains in its orbital and the final state of the 
He + ion is modeled as 

\^(E)) = 2-4 {\ls z 'ME)) + \4>c{E)ls z <)) (12) 

where \4> C (E)) is the Coulombic continuum wavefunction 
for energy E and charge Z = +1. Using these model 
wavefunctions the analytic expression 



P(Z e ,E) 



{l-Z e )Z° e 'e 



3/2 2 arctan(„/Z e ) 



v/l-e-^^ + Z 2 ) 2 



dE 



(13) 



with k = ^/{E + 2)/Ry and lRy = 13.60585972 eV is 
obtained for the probability density, i. e. for the model 
tail for the high-energy continuum states. A fit to the 
ab initio spectrum yielded Z e = 1.3074. As is evident 
from (|13p the probability density decays exponentially 
for high energies which is the most essential property. 
(In fact, it has been verified that all subsequent conclu- 
sions are unchanged, if a different model tail is used in 
which the initial-state charge is fixed to the one of the 
tritium nucleus (Z = 1) and the effective charge of the 
final Coulomb wave is used as fit parameter.) 

Fig. [2] compares the calculated final-state probability 
density with the analytical model tail and confirms the 
applicability of the latter for high energies, in fact already 
starting from about 200 to 250 eV, similarly as for T2. 
The integrated probability density (including the model 
tail for energies above 904 eV) yields a total probability 
of 29.65% for ionization of He following j3 decay of T — . 

The mean excitation energy E relative to the electronic 
ground-state energy of T ~~ can be obtained from 

E = J2 E n P n + 
n 

/■904eV poo 

I EP{E)dE+ / EP(E) dE. (14) 

^24.59cV J904cV 

Insertion of the final-state probability distribution cal- 
culated in this work in (TTJ]) yields E = 27.487eV for the 
mean excitation energy of the decay product. This result 
may be compared to the one obtained by the alternative 
relation 



(*r|H(He)|*r> 

£b(T-)-2<*ri-i*r>- 

r 



(15) 



G 




Energy [eV] 

FIG. 2: Final-state continuum probability density of He after 
/3 decay of T~ (solid) and a model tail with Z e = 1.3074 
(dashed) . The inset shows the probability density on a linear 
scale. 

With the expectation value of (r _1 ) and the ground state 
energy Eq(T~ ) reported by Frolov in [13] the mean ex- 
citation energy calculated with (|15[) is E = 27.469 eV. 
Frolov also reported the expectation values for the case 
of a finite mass of the tritium anion nuclei. Again us- 
ing dT5j) for this case one obtains E 1 = 27.479 eV. This 
comparison confirms the quality of the final-state distri- 
bution obtained in this work and validates furthermore 
the use of the approximation of an infinitely heavy nu- 
cleus approximation. 

Frolov noted in [l7| that his calculated bound-state 
probability appears to imply a continuum contribution of 
about 30%. Since this value is about 10 time larger than 
the continuum probability known for neutral T atoms, he 
speculated that in fact a large number of decays (about 
15 to 20 %) may end up in triplet states of helium, leaving 
a much smaller fraction in the singlet continuum. The 
reasoning in is that the more diffuse ground state 
of T - is not sufficient to explain such a large continuum 
contribution, as follows from a comparison to results ob- 
tained for Rydberg states of neutral T atoms. On the 
other hand, triplet states may be populated by the (vir- 
tual) interaction with the /3 electron omitted in the sud- 
den approximation. This argumentation is, however, er- 
roneous. As has been discussed in detail in [2l|, the sum 
rule for the sudden approximation gives always unity, in- 
dependently on higher-order corrections to it. This is 
also (as discussed above) fulfilled by the present calcu- 
lation which indeed confirms the about 30 % continuum 
probability indirectly found but rejected in [l7i ] . 

Finally, the exchange interaction is expected to be 
much smaller than the direct one (in [33] it was found 
for atomic tritium to be by a factor ij 2 smaller), but al- 



ready the direct term (first term beyond the sudden ap- 
proximation) is by a factor r/ 2 smaller than the sudden 
approximation. Close to the end point of tritium /? decay 
one finds for the Sommerfeld parameter 77 w —0.0271 and 
for T2 it was explicitly shown that the first-order correc- 
tion to the sudden approximation is itself of the order 
of r/ 2 and thus of the order of 0.01 % [llj . in accordance 
with corresponding system-independent sum rules given 

in [m. 

IV. CONCLUSION 

In this work the complete final-state probability distri- 
bution of He following the nuclear /3~ decay of tritium 
anions has been calculated. For the small number of 
bound states considered previously in [l7( the agreement 
is very good for the dominant ground and first excited 
states, but especially for the (very weakly populated) 
n = 3 state a deviation by about an order of magnitude 
was found. It is preliminarily attributed to a possible 
typographical error. The agreement with an earlier rel- 
ativistic multi-configuration Dirac-Fock calculation [lj[ 
is on the other hand very poor. Since such a large size 
of relativistic effects is not expected, especially not for 
the light nuclei involved, this deviation is attributed to 
a possibly too small basis set used in [lj[. Nevertheless, 
the present study may stimulate further theoretical work 
to clarify the discrepancies to [jj} and of both works to 
the relativistic ones in [l8j ]. 

In order to further test the accuracy of the present cal- 
culation, the bound-state energies were compared to very 
accurate literature data and were found to agree very ac- 
curately with them. Furthermore, the mean excitation 
energy obtained from the complete final-state spectrum 
was compared to the value predicted on the basis of clo- 
sure. Again, very good agreement was found. There- 
fore, the results of this work should be reliable and of 
direct importance for the tritium neutrino-mass experi- 
ment KATRIN that is presently under construction. In 
order to allow the use in the experimental analysis and 
for predicting how much a possible T~ admixture to the 
T2 source spoils the extracted neutrino mass, the contin- 
uum transition probability is given in binned form, but 
also available numerically on request. Finally, it may be 
noted that a controlled admixture of T~ to the tritium 
source may in fact be used for the analysis of the experi- 
mental sensitivity to the atomic and molecular final-state 
spectrum. 
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